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APPROXIMATE METHOD FOR CALCULATING FREE VIBRA- 


TIONS OF A LARGE -WIND-TURBINE TOWER STRUCTURE 
by Sankar C. Das* and Bradford S. Linscott 
Lewis Research Center 

SUMMARY 

An approximate method is presented to calculate the fundamental 
bending and torsion natural frequencies of a typical tower structure 
for a wind turbine. The tower structure is modeled as a system of 
six masses and six springs. This simplified lumped-mass model 
leads to a set of ordinary differential equations involving matrices in 
terms of mass, mass moments of inertia, structural stiffness, or 
structural flexibility. Solution of the differential equations leads to 
the common eigenvalue or characteristic value problem. 

Natural frequencies of the second and higher modes are often 
considerably greater than the fundamental natural frequencies for 
large frame structures. As a result, Dunkerley’s equation can be 
used as a solution to the eigenvalue problem and to determine the 
fundamental natural frequencies of the tower structure. 

The fundamental natural frequencies in bending and torsion were 
calculated for the MOD-0 wind turbine tower structure by the approxi- 
mate method. The frequencies were calculated for the tower with and 
without the nacelle and rotor blades. The approximate fundamental 
natural frequencies for the tower agree within 18 percent with test 
data and NASTRAN predictions. Successful design of wind turbine 
tower structures depends on the proper placement of the fundamental 
natural frequencies and modes. Thus, an approximate method for 

*Associate Professor of Civil Engineering, South Dakota State 
University, Brookings, South Dakota; Summer Faculty Fellow at the 
Lewis Research Center in 1976. 
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predicting tower frequencies should be used to aid in the preliminary 
design and sizing of tower members to meet frequency placement re- 
quirements. 

Once the tower design is finalized and the fundamental natural 
frequency placements determined with a detailed structural model, an 
approximate method can again be useful. It can be used to find gross 
errors in the detailed model, to differentiate between important modes 
and those caused by local member vibration, and to identify coupled 
mode shapes. Reasonable agreement between results obtained from 
the detailed model and those obtained from an approximate method 
allows the detailed model to be used with greater confidence in eval- 
uating the dynamic response of the coupled wind turbine generating 
system. This analysis is a key step in the program goal of predicting 
the dynamic response of wind turbines. 

INTRODUCTION 

Recent shortages in the supply of energy, coupled with increasing 
fuel costs, have forced our Nation to research all forms of energy. 

The power available from the wind is now being examined with renewed 
interest. 

The Federal wind energy program is directed by the Energy Re- 
search and Development Administration (ERDA). The program in- 
cludes research and development on a variety of applications and con- 
cepts for wind energy systems. Agreement was reached that, under 
the overall program management of ERDA, the NASA Lewis Research 
Center would provide project management for a portion of the overall 
program. 

As part of this program, Lewis has designed and constructed a 
wind turbine large enough to assess the technology requirements and 
the associated operational problems of a large wind turbine. The 100- 
kilowatt wind turbine has been constructed at the NASA Lewis Plum 
Brook Station near Sandusky, Ohio. The wind turbine consists of an 
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open-truss steel tower, 93 feet high; a nacelle that houses the alter- 
nator, the gearbox, and the low-speed drive shaft; and two aluminum 
rotor blades, each 62.5 feet long. The wind turbine is designed to 
produce 100 kilowatts of electric power in an 18-mph wind at a rotor 
speed of 40 rpm. 

For proper design of a wind turbine, it is necessary to perform 
analyses to determine the natural frequencies and associated mode 
shapes of the primary structures of a wind turbine, such as the tower. 
This report provides an approximate method for calculating the first 
cantilever -bending natural frequency and the first torsional natural 
frequency of a typical tower structure for a large wind turbine. 

A typical wind turbine tower structure can also be modeled and 
analyzed with the aid of computer programs such as STRUDL n or 
NASTRAN. However, as the tower model becomes more elaborate, 
the time required to prepare the program for running becomes lengthy. 
The purpose of this analysis is to provide an approximate modeling 
method that can be used in the preliminary design stage to quickly 
calculate the fundamental natural frequencies of a typical wind turbine 
tower. 

This report describes the mathematical model selected to re- 
present a typical wind turbine tower structure. The simplified model 
can be used to compute tower natural frequencies with the aid of a 
pocket-type calculator. The model can also be programmed with 
NASTRAN or STRUDL II to eliminate some of the need for repetitious 
calculations. The equations of motion for the tower structure are 
developed and solved in closed form. The equations for the free 
flexural vibrations of the tower are also presented. Finally, in an 
appendix, the natural frequencies are determined for the ERDA-NASA 
MOD-0 wind turbine tower structure by using the approximate method. 
The frequency values calculated by the approximate method are com- 
pared with values determined by testing. In addition, these funda- 
mental frequencies are compared with values obtained by analyzing a 
relatively much larger structural dynamic model of the tower with 
the aid of NASTRAN. 
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SYMBOLS 

cross-sectional area 

cross-sectional area of diagonal brace member 

flexibility coefficient in bending 
flexibility coefficient in torsion 
horizontal distance; derivative 
modulus of elasticity 
force 

frequency, co/27t 

shear modulus 

vertical distance 

area moment of inertia 

whole numbers 

polar mass moment of inertia 

stiffness coefficient 

length of diagonal brace member 

mass 

axial load 

generalized coordinate 
radial distance 

equivalent torsional rigidity of tubular tower legs and 
bracing (ref. 3) 
h/n 

kinetic energy 
time 

potential energy 
width 

Cartesian coordinates 
linear displacement 
cylindrical coordinate 
projected angle on y,z plane 
natural frequency 

first derivative with respect to time 
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(") second derivative with respect to time 

1, 2, 3, 4 tower sections 


Superscript: 

D 

0 

1 , 2 , 3 , 4 , 5, 6 


diagonal brace member 
center of tower 
tower bays 


ANALYSIS 

Mathematical Model Description 

The wind turbine tower structure is composed of four main col- 
umns with diagonal and horizontal bracing, as shown in figure 1. Be- 
cause the tower height is much greater than its average width, the 
tower is assumed to behave like a cantilever beam. 

The mathematical model for the tower treats the mass properties 
and the elastic properties separately. The tower mass is accounted 
for by lumping portions of the total mass at node points. A typical 
node point is located at the intersection of a horizontal member, a 
diagonal brace, and the vertical leg, as shown in figure 1. The four 
nodal masses at each tower level are added together, resulting in a 
single value for each of the six bays or levels of the tower. The cen- 
ter of each lumped mass is located on the vertical centerline of the 
tower and at the same elevation of each of the horizontal structural 
members. 

The elastic properties for the tower structure are calculated at 
various levels along the tower height. The elastic properties are 
determined by calculating the area moment of inertia of the structure 
at each tower level selected and by using the material mechanical pro- 
perties. 

The tower structure is thereby simplified to a system model com- 
posed of six lumped masses and six elastic beams. Figure 2 shows 
the tower structural model that is used to calculate the natural fre- 



quency of the tower in bending. The model used to calculate the 
natural frequency of the tower in torsion is shown in figure 3. 


Equations of Motion 

The equations of Lagrange (ref. 1) are used to formulate the 
equations of motion for the tower structural model. First, the kinetic 
energy of the multidegree system is expressed as 


n 


T = 


i=l 


1 M.(y )2 + 1 m.(zJ 2 + 1 j (0 )2 

2 ^ 2 ^^ 2 ^_ 


( 1 ) 


Next, the potential energy of the multidegree system, neglecting the 
effects due to gravity, is written as 


V 


n 

" / 
ti 



From the Lagrai^e equation. 


d /9T\ 3T ^ av 
dt l^aqj 3q. 3q. 


( 2 ) 


(3) 


where q^ represents the generalized coordinates y^, z^ and 6 ^ Equa- 
tions (1) to (3) yield 
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These equations are uncoupled linear ordinary differential equations of 
equilbrium in and 0^^. Therefore, free vibrations for the 

idealized structure in the coordinates y, z, and 6 are independent 
of each other. This is true when the mass center of the tower struc- 
ture coincides with the center of twist of the various tower levels . 

The system of linear differential equations (4) can be written in 
the matrix form 


[M] 

{y} 

+ [Ky] 

{yj =0 

(5) 

[M] 

{z} 


)z} = 0 

(6) 

[Jfl] 

{0} 


(0) =0 

(7) 


The mathematical model for the tower structure is symmetric for bend- 
ing in the y and z directions. As a result, the frequency solutions in 
the y and z directions will be identical. Since equations (5) and (6) 
yield identical results, only equation (6) is used in the following analysis. 

Free Flexural Vibrations' in z -Direction 

In equation (6), the term [K^] is the stiffness coefficient matrix for 
the structure in the z -coordinate direction. Often it is useful to gen- 
erate flexibility coefficients for the structure. Equation (6) may be 
rewritten, using the flexibility coefficient matrix [a^,] rather than the 
stiffness coefficient matrix [K^], as 

[a^] [M] {z}+{z}=0 (8) 




where 
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and 

[KJ = [I] 

the identify matrix. 

Equation (8) forms a set of linear second-order differential equa- 
tions whose solution is given by (ref. 1) 

{ a } = -w^ { z } (9) 

From equations (8) and (9), 

[M] {z}={z} (10) 


or, expanding, 



n 

Dividing equation (10) by w yields 



{z} = {0} 


( 11 ) 


( 12 ) 


Equation (12) forms a set of homogeneous algebraic equations in z. 
(1/w is unknown. ) The solution of these equations is usually called 
the eigenvalue or characteristics value problem. A nontrivial solution 
of equation (12) can exist only if the determinant 
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[al 


= 0 


( 13 ) 


w 



K 

^In " 


a, M„ 
" 


= 0 ( 14 ) 



nn 


where the mass matrix is a diagonal matrix, 



[M] = 






and the flexibility matrix in the z-coordinate direction is 
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^In 

21 

a . . 

^22 

• • • 3. 

‘^2n 






Assume the roots of the frequency equation (14) are l/w_ , 

^1 

l/w_ , 1/w^ , , l/w_ . Expanding equation (14) with the 

^2 ^3 ^n 

assumed roots leads to the following equation from reference 2 : 



. + 


1 




+ a_ M 
^22 


2 




Mo + . . 

33 


+ a„ 




nn 


n 


(15) 


Therefore, natural frequencies of the second and higher modes are 
often considerably greater than the fundamental frequency. If this 
holds, all terms on the left side of equation (15), except the first, may 
be omitted for the approximate determination of the fundamental fre- 
quency . It can be written as 








Mo 

33 


+ 


+ 


M_ 


nn 


n 


(16) 


Equation (16) is known as Dunkerley’s equation (ref. 2) and allows the 
fundamental frequency to be determined with reasonable accuracy by 
using longhand calculations. 
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Evaluation of Flexibility Coefficients 

By using the area moment principle suggested in reference 3, the 
flexibility coefficients can be written as 



k=n 



/ 

k=l 




(17) 


For the case of six degrees of freedom (i = 1, 2, . . . , 6 and j = 1, 

2 , . . . , 6 ), 



where h^ is the vertical distance between the tower node points and 
n is an arbitrary number of sections selected between node points . 

At each selected section, between node points, the area moment of 

Ir 

inertia I. is calculated for the tower about the y- or z-axis. 

^ It 

The bending moment Mj is calculated between two particular 
tower node-point levels (i = 1, 2, . . . , 6) and at a particular tower 
section (k = 1, 2, . . . , n). The magnitude of the bending moment 
is determined by multiplying a unit horizontal force applied at a 
node point by the vertical distance Xj^. The quantity x^ is the dis- 
tance between a section of the tower (k = 1, 2, . . . , n) where the 

tr 

bending moment intensity occurs and the unit force applied at a 
node point (j = 1, 2, . . . , 6). The quantity E is Young’s modulus 
of elasticity for the particular tower structural material. 

In appendix A of this report the first flexural frequency for the 
MOD-0 wind turbine tower structure is calculated by using the approxi- 
mate method as described herein. Each bay of the tower having a 
height hj was divided into n = 4 sections. At each of 24 sections 
(n = 4 sections times j = 6 bays) of the tower, the area moment of 
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inertia was calculated. The bending moment m| due to a unit hori- 
zontal force sequentially applied to each of the six node points, was 
computed at each section. Then equation (17) was used to calculate 
the flexibility coefficients a„ , where i = j for each of the six tower 

/ 

bays. Finally, equation (16) was used to calculate the first natural 
bending frequency for the MOD-0 wind turbine tower structure with 
and without the nacelle and rotor blades. 


Free Torsional Vibrations 

In equation (7) the term [K^] is the torsional stiffness coefficient 
matrix for the tower structure. Often it is useful to generate torsional 
flexibility coefficients for the structure. Equations (7) may be rewritten, 
by using the torsional flexibility coefficient matrix [b^] rather than the 
stiffness matrix [K^], as 

[b«] [J] {9}+{e}=0 (18) 

where 

[b^] [Kj,] = [I] (19) 


is the identify matrix. 

The solution of equation (18) is similar to that of equation (8) . 
Therefore, the fundamental torsional frequency can be written, by using 
Dunkerley’s equation, as 



b, 


'11 


+ b 





( 20 ) 


Evaluation of Torsional Flexibility Coefficient 

The torsional spring constant for the tower structure can be written 
as 
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8(Ajj) E 


^9 = 




1 2 2 
— Wj cos sin + 4 



^GR 


+ 4 


eqA 


“i 


( 21 ) 


where 

cross-sectional area of a diagonal brace member 

E modulus of elasticity 

vertical tower height between node points 

Wj half of width of a square horizontal frame 

angle of diagonal brace member projected on z-y plane 

Ij area moment of inertia of a tubular leg taken about centroid 

of tube 

R- radial distance from tower x-axis to center of a tubular leg 

G shear modulus 

RgQ equivalent torsional rigidity of tubular tower legs and bracing 
^ (ref. 3) 

The total torsional spring constant is composed of three terms 
in equation (2 1) . The first term accounts for the change in length of the 
diagonal crossmember bracing as located on the sides of the tower. 

A torsional moment vector applied along the y-axis of the tower re- 
sults in tension and compression loads in these crossmembers. (This 
term is developed in appendix B.) The second term in equation (21) 
accounts for beamwise bending of each of the four main tubular legs of 
the tower. Again, a torsional moment vector applied along the x-axis 
of the tower tends to rotate the square frame structure within a hori- 
zontal y-z plane. Rotation of each frame causes bending in the tubular 
tower legs. 

Rotation of each frame structure in the y-z plane also causes 
torsion in each tubular tower leg and in the diagonal cross bracing. 

The third term in equation (21) accounts for the torsional stiffness of 
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these members. For a typical wind turbine tower structure, such as 
the MOD-0, the magnitudes of the last two terms on the right in equa- 
tion (21) were found to be small as compared with the first term. As 
a result, the last two terms were neglected when the torsional stiff- 
nesses for the MOD-0 tower structure were calculated in appendix C. 

The first torsional frequency for the MOD-0 wind turbine tower 
structure was also calculated in appendix C, by the approximate method 
described herein. Equation (2 1) was used to compute the torsional 
stiffness coefficient at each of the six tower bays. The polar mass 
moment of inertia J was calculated at six tower levels. Finally, equa- 
tion (20) was used to calculate the first torsional frequency of the tower 
structure with and without the nacelle and rotor blades, also in appen- 
dix C. 


SUMMARY OF RESULTS 

Calculated and experimental MOD-0 tower frequencies are com- 
pared in tables I and II. Table I presents a summary of MOD-0 tower 
natural frequencies without the nacelle and rotor blades. The test data 
were taken from reference 4. The values calculated by using NASTRAN 
were taken from references 4 and 5. The approximate values in table I 
were calculated by using the approximate method described in this re- 
port. The detailed calculations are presented in appendixes A and C. 
The 5. 04-hertz first bending frequency calculated by the approximate 
method is 7 percent higher than the 4. 7-hertz value determined by test- 
ing. The 8. 59 -hertz first torisonal frequency calculated by the approxi- 
mate method is 18 percent lower than the 10. 5 -hertz value determined 
by testing. 

Table II summarizes MOD-0 tower natural frequencies with the 
nacelle and rotor blades mounted on top of the tower. The rotor blades 
were oriented in the horizontal position, parallel to the ground. The 
test data were taken from reference 4. The values predicted by using 
NASTRAN were taken from references 4 and 5. The approximate values 
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in table EE were calculated by using the approximate method described 
in this report. The detailed calculations are contained in appendixes A 
and C. The 2. 45-hertz first bending frequency calculated by the approxi- 
mate method is 15 percent higher than the 2. 1-hertz value determined 
by testing. The 3. 78 -hertz first torsional frequency calculated by the 
approximate method is 16 percent lower than the 4. 4-hertz value deter- 
mined by testing. 

As a result of this analysis, an error was found in references 4 
and 5. The first torsional tower frequency with the nacelle and rotor 
blades (horizontal) was reported as 9.8-hertz in reference 4 and 9.56- 
hertz in reference 5. This approximate analysis predicts the first 
torsional natural frequency as 3.78-hertz. Reference 4 reports that 
a 4. 4-hertz mode "appears” to be torsion, with first-mode bending 
present in the north-south direction. This approximate analysis veri- 
fies that the 4. 4-hertz mode is, in fact, the first torsional mode and 
that the 9. 8 -hertz mode is clearly not the first torsional mode. 

For structural models that have a large number of degrees of 
freedom, it is very difficult to determine certain fundamental modes 
and frequencies. This difficulty was experienced during the work re- 
ported in reference 5, since the first torsional tower frequency was 
established at 9. 5 6 -hertz. If an approximate frequency analysis, as 
presented herein, was first conducted, the fundamental modes and 
frequencies as predicted by NASTRAN could have been more easily 
and accurately interpreted. 

CONCLUSIONS 

The results of this analysis show that there are important reasons 
for using an approximate method to evaluate the fundamental natural 
frequencies of a typical wind turbine tower structure. The reasons 
are summarized as follows: 
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1. The simplified model can be used during preliminary tower 
design to easily determine tower fundamental frequency placement. 

2 . The simplified model can be used to check frequency predic- 
tions resulting from the analysis of a more-detailed structural model 
of the tower. A detailed structural dynamic model is needed to assist 
in finalizing the tower design. 

These checks include 

a. Distinguishing between overall tower system frequencies 

and individual tower member frequencies 

b. Identifying coupled modes 

c. Cheekily for potential gross frequency prediction errors 
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APPENDIX A 

APPROXIMATE CALCULATION OF FIRST FLEXURAL 
NATURAL FREQUENCY FOR MOD-0 TOWER 
Procedure for Computing Flexibility Coefficient 

1 of Tower 

The procedure for computing the flexibility coefficient for 
tower bay 1 is as follows: Rewriting equation (17) from the ANALYSIS 
gives 




h-1 


.1 


Specializing this equation for tower bay „l, shown in figure <1, yields 


'll 


h / . AA 


E 


3 3 4 

mJxJ M]xJ 


I 




L 


First, computing the values for the bending moments at each level 
shown in figure 4 gives the bending moment due to the horizontal 1- 
pound force about level 1 as 

mJ = 1.0 Ibx 31.5 in. = 31.5 in-lb/lb 


then 
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= 94. 5 in-lb/lb 

=157.5 in-lb/lb 
Mj = 220. 5 in-lb/lb 

The distances from the point where deflection due to unit load are 
desired to the locations of the bending moment intensity are 

xj = 31. 5 in. 

Xj = 94. 5 in. 

Xj = 157. 5 in. 

X2 = 220. 5 in. 

Each of the four vertical legs of the tower (in bays 1 and 2) is fabri- 
cated from 8. 0 -inch-diameter extra-heavy pipe. The cross-sectional 

area of the pipe is 12. 76 square inches; the area moment of inertia of 

4 

the pipe is 105. 7 inches. 

The transfer formula for computing the area moment of inertia of 
the pipe about the tower center is 

I = Iq + Ad^ 


where 


Iq = 105. 7 in'^ 


A = 12.76 in^ 


and 
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and d is the horizontal distance from the center of the tower to the 
center of the tubular tower leg. 

Then, from figure 4, accounting for four tubular legs, the area 
moments of inertia for bay 1, sections 4 to 1, are 

= 4 [l05. 7 + 12. 76(134)^ = 9. 1758X10^ in^ 

= 4 [l05. 7 + 12. 76(147. 2)^ = 1. 1060X10 ^.in'^ 

= 4 [^5. 7 + 12. 76(160. 3)^ = 1. 3121X10® in'^ 

I J = 4 [io5. 7 + 12 . 76(173. 4)^ = 1. 5358X10® in'^ 

Again referring to figure 4, Sj = 64 in. 

Finally, substituting the values for Mp x, and I into the 
equation for the flexibility coefficient yields 

„ _ 63r(31.5)(31.5) 

^11 T 

^ |_9. 1758X10^ 

= MM in/lb 
■ E 

where E is the material modulus of elasticity. 

Procedure for Computing Flexibility Coefficient ^22 
for Bays 1 and 2 of Tower 

The procedure for computing the flexibility coefficient a .22 for 
tower bays 1 and 2 is as follows: 


(94. 5)(94.5) ^ (157.5)(157.5) ^ (220. 5)(220. 5) 
1.106X10® 1.3121X10® 1.5358X10® 
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^22 " 


E / Tk E 


"^2 2 


k=l 


k=l 


where 


= 25. 5 + (k - 1)S, k = 1, . . . , 4 


Mj = = 204 + j Sj + (k - l)Sj k = 1, . . . , 4 


The values for I, and S are given in table in. Subsituting these 
values into the preceding equation yields 


a22 =M:^ln/lb 
E 


Procedure for Computing Flexibility Coefficient a^g 
for Bays 1 to 3 of Tower 


The procedure for computing the flexibility coefficient agg for 
tower bays 1 to 3 is as follows: 


^33 




4 



! 

k=l 



where 
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Mg = = 24 + (k - 1)S^ ■ k = 1, . . . , 4 


= 192 + -^ + (k - 1)82 k = 1, . . . , 4 
E 


Mj = 3^ = 396 + — + (k - 1)8 j k = 1, . . . , 4 

E 

The values for I and 8 are given in table ni. 8ubstituting these 
values into the preceding equation yields 


a33 = il^ln/lb 


E 


Procedure for Computing Flexibility Coefficient a^^ 
for Bays 1 to 4 of Tower 

The procedure for calculating the flexibility coefficient a^^ for 
tower bays 1 to 4 is as follows: 


4 4 

8.\ ' mW 8, 


^44 " 


Z__. 

k=l ^ 


E 


,k E 


Tk E 


k=l 


k=l 


k=l 


M^x^ 


I 


M4 = X4 = 21 + (k - 1)84 k = 1, . . . , 4 

S 

= x^ = 168 +--^ + (k - 1)8« k = 1, . . . , 4 


where 
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Mg = Xg = 324 + ^ + (k - l)Sg 

g 

Mg = x^ = 516 + -2 -H- :(k - IjSg 
2 

S 

= x^ = 720 + + (k - 1)S^ 

2 


k= 1, . . . , 4 
k = 1, . . . , 4 
k = 1, . . , 4 


The values for I and S are given in tables in and IV. Substituting 
these values into the preceding equation yields 

agg = in/lb 

00 g 


Procedure for Computing Flexibility Coefficient a^g 
for Bays 1 to 6 of Tower 

The flexibility coefficient agg for tower bays 1 to 6 is calculated 
as follows: 




'66 


E 


k=l 


E 


k=l 


^ i; 


+ 

k E 


M^St^ 

ivi3Xg 


k=l 


_/ r 


S^X Sg 

T / + — 


4_^ 4 




mV 
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k=l 


i! 
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k=l 


y I, 


k E 
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k=l 


where 
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+ 

00 

II 

II 

T 

- l)Sg k 

II 

, 4 

= 144 + — + 

(k - DSj 

k= 1 , . 

. . , 4 

Ic It ^4 

= 300 + — + 
2 

(k - DS 4 

k= 1 , . 

. . , 4 

Mg = x^ = 468 + — + 
2 

- DSj 

k= 1 , . 

. . , 4 

It Ic ^9 

= X 2 = 660 + + 

2 

Oc - DSj 

k= 1 , . 

. . , 4 

Mj = x^ = 864 + — + 

(k - l)Sj 

k= 1 , . 

. . , 4 


The values for I and S are given in tables in and IV. Substituting 
these values into the preceding equation yields 


= Mil in/lb 


"66 


E 


The tower structure weight is estimated at 44 000 pounds. The total 
tower height is 93 feet. Multiplying the ratio of bay height to total 
tower height times the total mass of the tower structure yields the 
estimated mass of each tower bay. Then the mass for bay 1 is 


M. = ^ = 25. 71 lb-sec^ /in 

^ (32.2)(12) 93 


where the height of bay 1 is 21 feet, as shown in figure 4. For the re- 
maining tower bays the masses are 



25 


M 2 = 113.87^-^^= 20.82 Ib-sec-Vin 
M 3 = 113. 87 = 19. 59 Ib-sec^/in 

= 113. 87 = 17. 14 Ib-secVin 

Mg = 113. 87 ^ = 15. 92 Ib-sec^/in 

Mg = 113. 87 = 14. 70 Ib-sec^/in 


Rewriting Dunkerley’s equation gives 
1 




^ 11^1 ^ 22^2 ^ 33^3 ^ 44^4 ^ 55^5 ^ 66^6 


Substituting the values for the masses and the flexibility coefficients 
then yields 


1 ^ 1 


(^i) 


2 E L 


(3. 76)(25. 71) + (26. 64)(20. 82) + (97. 36)(19. 59) 


+ (252.)(17. 14) + (535. 7)(15. 92) + (986. 2)(14. 




29 904 


.6 


where E is 30x10 psi. The frequency is therefore estimated as 

f 6 \ 1/2 

^ 277 ^ 277 \29 904 / 
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The tower first bending frequency is reduced by the addition of 
the stairway, the nacelle, and the rotor blades. These components 
of the wind turbine are shown in reference 4, with their associated 
weights. The tower frequency with these components added is easily 
estimated. The stair weight (12 000 Ibm) is divided into six equal 
weights and lumped at each node point. The nacelle and rotor blade 
weight (34 000 Ibm) is lumped at node 6. Then the new masses are 

M. = 25.7+ 25.7 + 5.2 = 30.9 lb-sec^ /in 

^ (32.2)(12)V6/ 

M2 = 20. 8 + 5. 2 = 26. 0 Ib-sec^/in 

Mg = 19. 6 + 5. 2 = 24. 8 Ib-sec^/in 

M^ = 17. 1 + 5.2 = 22. 3 lb-sec^ /in 

Mg = 15.9 + 5.2 = 21. 1 Ib-sec^/in 

Finally, 

M(.= 14. 7 + 5. 2 + = 19. 9 + 88 = 107. 9 lb-sec^ /in 

° (12)(32.2) 


Substituting the new values for the masses and the same values 
for the flexibility coefficients into Dunkerley’ s equation yields 


1 


]_ 

E 


(3. 76)(30. 9) + (26. 64)(26. ) + (97. 36)(24. 8) 


+ (252)(22. 3) + (535. 7)(21. 1) + (986. 2)(107. 9) 



E 


126 600 
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The frequency is therefore 


2tt 2ti \126 600/ 


2.45 Hz 
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APPENDIX B 


DEVELOPMENT OF TORSIONAL SPRING CONSTANT 
FOR MOD-0 TOWER 


The derive the torsional spring constant for a typical bay of 
the tower structure, four horizontal forces F are applied to the node 
points of the tower structure as shown in figure 5. This set of forces 
cause a rotation 6 in the horizontal plane . Rotation is allowed by 
extension of the diagonal members a^bg, ^2^3’ ^3^2’ 

Extension of a single diagonal member may be expressed as 



^D^D 


^D^D 


where 

^D axial load in diagonal member 

length of diagonal member 

^D cross-sectional area of diagonal member 

Ej-j material modulus of elasticity for diagonal member 

From figure 5, 6^^ = R0, Sg = R6/\/2, and w = R/ \/2. Then 5^ = w0. 

As shown in figure 6 the horizontal force component of the axial 
load Pj-j in the diagonal member is expressed as F = Pj^ cos <p or 


^d-^d^d 

F = cos (p 


I 


D 


Again from figure 6, = §2 cos q> =vi cos ip. Then 



From figure 5, Zj) =h/sin <p. Then 


2 

F = w0 cos (p sin (p 

h 


For unit values of F and 9, 


„ 2 
F = w cos (p sin (p 


Next, accounting for the existence of two diagonals on one side, one 
carrying a load in tension and the other carrying an equal but opposite 
load in compression. 


2AdEd 2 

F = w cos cp sin cp 

h 


Finally, the unit value of the torsional moment M is written as 
M = 4(Fw) or 


2 2 

M = w^ cos^^ sin cp 


The unit value of M is also the torsional spring constant K^. 
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APPENDIX C 

APPROXIMATE CALCULATION OF FIRST TORSIONAL 

NATURAL FREQUENCY FOR MOD-0 TOWER 

From the equation for developed in appendix B and the struc- 
tural properties listed in table V,. the values for , . . . , 

Kn are calculated. Then the torsional flexibility coefficients are 
^6 


bii = 


13.5X10^ 


= 7. 407X10“^^ rad/in-lb 


b 22 = — + — = 22. 96X10"^^ rad/in-lb 


^1 h 


bgg = + J_ = 32. 39X10"^^ rad/in-lb 


K. K. K. 
^1 ^2 ^3 


b. . = — + — + — + — = 45. 24X10"^^ rad/in-lb 
K. K. K. Kn 

h h h ^4 

b.p. = — -t- — + — + — + — = 65. 48X10"^^ rad/in-lb 
K. K. K. K. K. 

^1 ^2 ^3 ^4 ^5 


and 


b(,„ = -1— + _L_ + -JL. + = 97 . 02x10"^^ rad/in-lb 

Kn Kn Kn Kn Kn Kn 

^1 ^2 ®3 ^4 ^5 ^6 
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Calculate Polar Mass Moments of Inertia 

Next the lumped polar mass moments of inertia Jg, . . . , J0 
must be calculated for each bay of the tower. The lumped masses 
Mj, Mg, • • • , Mg for each bay of the tower had been previously 
calculated. The distances d^, dg, • • • , dg from the center of the 
tower to a node point are shown in figure 7. 

Each mass is assumed to be lumped at a node point on the tower. 
Then 


= 2Mj(dj)^ = 2(25.7)(127. 5)^ = 8. 356x10^ Ib-in-sec^ 

Jg = 2Mg(d2)^ = 2(20.8)(85)^ = 3.006x10^ Ib-in-sec^ 

Jg = 2Mg(d3)2 = 2(19. 6)(69. 5)^ = 1. 893x10^ Ib-in-sec^ 
2M^(d'^)2 = 2(17. 1)(59;)2:= 1. 190X10^ Ib-in-sec^ 

Jg = 2M5;(dg)2 = 2(15. 9)(49.25)2 = 7. 713X10^ Ib-in-sec^ 

Jg = 2Mg(dg)^ = 2(14. 7) (40. 25)2 ^ 4.763X10^ Ib-in-sec^ 

Calculate Polar Mass Moment of Inertia for Nacelle Mounted on 

Tower with Rotor Blades Oriented Vertically and Horizontally 

The polar mass moment of inertia with the blades mounted verti- 
cally (fig. 8) is 
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Md^ = __ 3Q ..Q 0 . 0 ... , (75)2 ^ — 4000_ ^gg)2 
32.2 X 12 32.2 X 12 


= 436 724 + 281 832 = 718 556 Ib-in-sec^ 


With the blades mounted horizontally (fig. 8), 


Jh = 



Md^ 


(75)2 ^ — 4000 — (105)2 ^ — 4000 — (270)^ 
32.2 X 12 32.2 X 12 32.2 X 12 


= 1 437 214 in-lb-sec" 


Rewriting Dunkerley’ s equation gives 


(^1)^ 


~ ^ll'^l ^22*^2 ^33'^3 ^44*^4 ^55*^5 ’^66‘^6 


Substituting the values for the torsional flexibility coefficients and the 
polar mass moments of inertia then yields 

— i — = 7. 407X10"^^ (8. 356X10^) + 22. 96X10"^^ (3. 006X10^) 


+ 32.39X10"^^ (1.893X10^) + 45.24x10"^^ (1. 19x10^) 
+ 65. 48X10"^^ (7. 713X10^) + 97. 02X10"^^ (4. 763X10'^) 
= 3:4278X10"^ 

(c0j)2 s 0.29175X10^ 

= 54. 0 rad/sec 
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The first torsional frequency is therefore 

f = -L CO, ^ Mil = 8. 59 Hz 
277 ^ 27T 

Calculate First Torsional Frequency of Tower with 

Nacelle and Blades Mounted Vertically 

Modifying the last term of Dunkerley’s equation to account for the 
polar mass moment of inertia of the nacelle and blades in the vertical 
position gives 

bgg(Jg + Jy) = 97. 02X10"^^ (4. 763X10^ + 7. 185x10^) 

Then 

— = 3 . 42 78x10“^ + bggJy = 3. 4278x10“'^ + 97. 02x10"^^ (7. 185x10^) 
= 10. 399x10"'^ 

= 9.616X10^ 
w = 31. 0 rad/sec 

The first torsional frequency is therefore 

f = J_ (co) = ii = 4. 93 Hz 
2t7 2tt 
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Calculate First Torsional Frequency of Tower with 

Nacelle and Blades Mounted Horizontally 

Modifying the last term of Dunkerley’s equation to account for the 
polar mass moment of inertia of the nacelle and blades in the horizontal 
position gives 

bgg(Jg + Jjj) = 97. 02X10"^^ (4. 763X10^ + 1. 4732X10®) 

Then 

— — = 3. 4278X10“^ + bggJjj = 3. 4278X10"'^ + 97. 02X10"^^ (1. 4732X10®) 

(C0j)2 

= 17.7208X10"^ 

(c0j)2 = 5.6431 x102 
= 23. 75 rad/sec 

The first torsional frequency is therefore 

f = = 23^ ^ 3 

2ir 


2ir 
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TABLE I. - SUMMARY OF MOD-0 TOWER NATURAL FREQUENCIES 
WITHOUT NACELLE AND ROTOR BLADES 


Tower mode 

Approximate 

values 

”1 

Test 

data^ 

NASTRAN 

predictions^ 

Tower frequency, Hz 

First bending in north -south direction 
First torsional 

5.04 

8.59 


4 . 76 
10 1 


^From reference 4. 


TABLE II. - SUMMARY OF MOD-0 TOWER NATURAL FREQUENCIES 
WITH NACELLE AND ROTOR BLADES ON TOWER IN 
HORIZONTAL POSITION 


Tower mode 

Approximate 

Test 

NASTRAN 


values 



data^ 

predictions^ 


Tower frequency, Hz 

First bending in north-south direction 

2.45 

2. 1 

2- 15 

First torsional 

3.78 

4 4^ 

(c) 


^From reference 4. 

^Erronenously reported as 9. 8 Hz in reference 4. 
^Erronenously reported as 9. 56 Hz in references 4 and 5- 
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Bay Sectiqn Area moment Cross-sectional Horizontal Area moment S 


of inertia 
at center 
of tower, 

lo’ 


area of center 
of tower, 
Aq, 


distance, 

d, 

in. 


of inertia, 

I, 



134.0 

9. 1758X10 

147.2 

1. 1060X10 

160. 3 

1. 3121 

173.4 

1. 5358 

122.2 

7.6244x10' 

111.6 

6. 3568 

100. 9 

5.2044 

90. 3 

4. 1672 

83.062 

2. 3208X10' 

79. 187 

2. 1096 

75. 312 

1. 9084 

71.437 

1. 7174 
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TABLE IV. - MOD-0 TOWER STRUCTURE SECTION 
PROPERTIES - BAYS 4 TO 6 


Bay 

Section 

Area moment 
of inertia 
at center 
of tower, 

io> 

in^ 

Cross -sectional 

area of center 

of tower, 

Aq, 

. 2 
in 

1 

Horizontal 

distance, 

d, 

in. 

Area moment 
of inertia, 

I, 

in^ 

S 

4 

4 

72 

.5 

8. 

40 

68. 187 

1. 5649x10^ 

42 


3 





65.562 

1. 4470 




2 





62.937 

1. 3337 




1 




60. 312 

1.2250 



5 

4 ' 

72 

.5 

8. 

40 

57. 781 

1. 1246X10^ 

39 


3 





55. 349 

L0321 




2 





52.906 

9.4327X10'^ 




1 


f 


1 

50.469 

8.5862 



6 

4 

I 72 

.5 

8. 

40 

i 48. 125 

7. 8099X10^ 

1 36 


3 

j 




1 45.875 

7.0993 




2 





! 43.625 

6.4228 




1 


' 



41. 375 

1 

5. 7803 
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TABLE V. - DIMENSIONS AND TORSIONAL SPRING CONSTANTS 


FOR EACH TOWER BAY 


[Modulus of elasticity, E, 30x10^ Ib/in^. ] 


Bay 

Cross-sectional 
area at center 
of tower, 
Ad? 
m2 

Width, 

w, 

in. 

Vertical 

distance, 

h, 

in. 

' ■ 1 ■ . 
Projected 
angle on 
y,z plane, 

<P, 

deg 

Spring constants, ^ 
K^, 

in-lb/deg 

1 

3. 17 

127. 5 

252 

54.5 

13.50X10^ 

2 

3. 17 

85. 

204 

58.0 

6.43;' 

3 

5.72 

69.5 

192 

51.2 

10.60 

4 

5. 34 

59.0 

168 

52.6 

7.78 

5 

4. 96 

49 ..2 

156 

55.2 

4.94 

6 

4.96 

40.2 

144 

58. 1 

3. 17 


% 


2 2 

wf cos (p^ sin cp^, where i = 1, 2, 
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Figure 2. - Model of tower 
structure for flexural 
vibration - six degrees 
of freedom. 
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Figure 3 - Model of tower 
structure for torsional 
vibration - six degrees 
of freedom. 










Figure 6. - Displacement of typical diagonal Figure 7. - Arrangement of tower 

brace member at node point. structure. (Dimensions are in 

inches. I 


Centerline 
of blade 


I 

Center 
of towers 

\ 

V 


I 


Nacelle weight 
without blades, 
30 000 Ibiri7 


I 



I 



■132 ► \ 


270 


Centerline 
of nacelle 


33— H 


-Blade weight. 
2000 Ibm 
-Centerline 
of blade 


Figure 8. - Plan view of centers of gravity and locations of MOD-0 
wind turbine components. (Dimensions are in inches. ) 
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